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1. Introduction 

1.1. Vacuum energy density and cylinder kernels 

The energy density of a quantum field in its vacuum state (or a state of fixed 
temperature), given particular time- independent external conditions, is of interest in 
cosmology, hadronic physics, and soft condensed matter physics (where the fluctuations 
involved are of thermal rather than quantum origin). Most famously, electromagnetic 
vacuum energy gives rise to the observed Casimir force between conducting bodies (and 
its counterpart for dielectrics, the Lifshitz force) Jj] The study of gravitational effects 
(or their negligibility, even in the presence of divergences) requires detailed knowledge 
of all components of the stress-energy-momentum tensor of the field as functions of 
position in space. This local information is also helpful in addressing such subjects 
of continued investigation and debate as the occasionally counterintuitive signs of the 
calculated forces, and the meaning of the divergences resulting from idealized boundary 
conditions. 

A convenient way to study vacuum energy is to insert (usually temporarily) an 
exponential ultraviolet cutoff in the spectral expansions of the stress tensor or the total 
energy. For the fiducial toy model, a scalar field, this procedure leads quickly to the 
study of the cylinder (or Poisson) kernels, which are Green functions for certain elliptic 
boundary value problems in one higher spatial dimension [231 [5j [10]. In cavities or 
billiards the cylinder kernel lends itself to the same kind of semiclassical or optical 
expansion that is more often (e.g., [29]) applied to the Green functions for the heat, 
Schrodinger, or resolvent kernels. (I hope to publish elsewhere a detailed justification 
and development of this assertion.) For rectangular cavities this construction reduces to 
the classic "method of images" and yields exact results [TlJ [21] , revealing, for example, 
that the divergent pressure on one face of a rectangle is related to the divergent energy 
density parallel to the adjacent faces. This insight is clearly pertinent to the vexed 
question of the vacuum pressure on a conducting sphere, where perpendicular boundary 
surfaces do not exist. Unfortunately, in the presence of curved surfaces (or, a fortiori, 
edges and corners) the optical construction is no longer exact. Preliminary calculations 
show that the leading semiclassical approximation is not adequate to understand the 
radial pressure in a sphere, and that in problems with curved boundaries the construction 
of higher-order approximations is blocked by problems of the same type that arose in 
connection with the heat kernel [251 El]- When the boundary is not smooth, the leading 
optical approximation is completely unacceptable until supplemented by diffractive 
contributions (e.g., [TF[ [7]), which are hard to calculate in generality. 

\ In lieu of a long list of references on vacuum energy, I recommend the recent special issues of New 
Journal of Physics [33] and Journal of Physics A [34j and the bibliographies of the papers therein. 
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1.2. Boundary integrals and multiple reflection 

Although higher-order optical approximations for smooth curved boundaries and 
diffractive approximations for nonsmooth ones can sometimes be found by trial and error 
and verified a posteriori, a more systematic approach is desirable. In principle, such an 
approach is available in the standard mathematical literature on partial differential and 
integral equations [22j [301 EH El EH]. A partial differential equation with boundary 
conditions can be reduced to an integral equation on the boundary. Up to some 



convergence issues touched upon lightly in Appendix B , the solution of the integral 
equation by the methods of Volterra or Fredholm is constructive. The result is a multiple 
reflection expansion (MRE) expressing the Green function of the original problem as 
a series J2n=qGn , wnere Gn is an iV-fold integral over the boundary; the solution 
is formally a "sum" over all paths from the source point to the field point, each path 
bouncing off the boundary iV times (in general nonspecularly). This construction was 
famously applied to physical problems by Balian and Bloch [H [2, Ej. (See also [26] 
in the case of the heat equation.) The more familiar (and simpler) semiclassical or 
optical approximations (involving only specular reflections) for the heat, Schrodinger, 
and resolvent equations emerge from the MRE when the integrals are approximated by 
steepest descent or stationary phase. In the case of the cylinder kernel treated directly 
by MRE, the appropriate approximation method is not so obvious; cylinder kernels can 
also be obtained from the other kernels by integral transformations, but at the cost 
of integrating parameters over values where the validity of the (e.g.) stationary-phase 
approximation is dubious. This circumstance calls for a careful examination of the 
implications of the MRE for the representation and approximation of cylinder kernels 
for curved or nonsmooth boundaries — and, indeed, at the first step, even for flat smooth 
boundaries. 



1.3. Quantum graphs 

The present paper is a preliminary foray with primarily pedagogical intent. The 
application of the MRE to vacuum energy needs to be studied first in the simplest 
case, one spatial dimension. But a one-dimensional cavity is just an interval, too 
trivial to occupy one's attention for long. A generalization, which combines the ease of 
exact solution of a one-dimensional problem with some of the nontrivial properties of 
multidimensional systems, is the concept of a quantum graph. 

A quantum graph J20J[T5] is a Riemannian one-complex — that is, a network of edges 
and vertices equipped with a natural notion of arc length on each edge. There follows 
a natural definition of the Laplacian operator as d 2 jdx 2 on each edge, supplemented by 
boundary conditions at each vertex to make the operator self-adjoint. In particular, the 
most natural and popular boundary conditions are the Kirchhoff conditions: Functions 
in the domain of the operator are required to be continuous at the vertices, while the 
outward derivatives of a function on all the edges incident on a given vertex are required 
to sum to zero. (A vertex attached to only one edge is thus a Neumann endpoint in the 
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usual sense.) 

On a Kirchhoff quantum graph the semiclassical approximation again reduces to the 
(exact) method of images, which is, nevertheless, quite intricate to execute. Solutions 
for the heat and resolvent equations have been known for some time [271 [T9l HI [18]. The 
corresponding construction for the cylinder kernel has been extensively investigated by 
Wilson et al. [32], [TJl [13l [HI [6]. Here I show how the same results can be obtained from 
an MRE, temporarily ignoring the fact that the simpler ansatz is already exact. The 
connection is not exactly trivial, so the exercise will be useful in tackling more serious 
models later. (It also enables us to get out of the way some complications peculiar to the 
one-dimensional case, associated with the logarithmic form of the integral kernel (CQ).) 

Section [2] sets up some machinery. The next four sections consider four increasingly 
complicated scenarios: an interval with one endpoint, a graph with one vertex but 
arbitrarily many edges, an interval with two endpoints, and finally a general quantum 
graph (with finitely many edges, finite in length). In each case the MRE calculations 
are carried only to the point where it is obvious how to continue them and match them 
up with known results. The main point becomes visible in the third case (section [5]), 
where a nontrivial boundary integral equation must be derived and iteratively solved 



for the first time. Some necessary integrals are evaluated in Appendix A Of greater 



interest is Appendix B where the mathematical status of the multiple-reflection series 



is discussed, and the possible utility of the exact Fredholm solution in the situation of 
finite temperature is pointed out. The machinery for finite temperature is set up in 



Appendix C 



2. Notation 

Let T be a quantum graph, and let V = R x T be the corresponding "Euclideanized 
space-time". Variables x, y, ... stand for points in T and variables s, t, ... for real 
numbers. I consider here only the standard Kirchhoff boundary conditions. 
The Green function for the Poisson equation in R 2 is 

1 ( r 

G (t,x;s,y) = --ln{- 

= -^H{t-sf + {x~yf}+C. (1) 

Let G(t, x; s, y) be the Green function for the Poisson equation in T. Then the cylinder 
kernels in T are 

T{t,x,y) = -2G{t,x;0,y), (2) 
dG 

T{t 1 x,y) = -2—(t,x;0,y). (3) 

Go is defined only modulo the indicated scale ambiguity (the arbitrary constants r and 
C), but Tq and the spatial derivatives of Gq are unique. Henceforth we take C — 0. 
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For later use note 

***( t x-s v) = - — — = (4) 

dx { , , ,V) 2n (t — s) 2 + (x — y) 2 dy ' { ) 

Related to this function is the well known distributional identity 
1 x 

lim— — - = 5(t — s), (5) 

xlO TT (t- S) 2 +X 2 V ' K ' 

which will be used repeatedly in what follows. 
3. Case 1: T = R + 

We will solve the Poisson equation on the half-plane with Neumann boundary at y = 0. 
Make the ansatz 

G = G Q + 7 , (6) 

/oo 
dsG (t,x;s, 0)/i (s), (7) 
-oo 

where /io has a dependence on (s , Vo) that will be notationally suppressed. The integral 
in ([7]) is over the boundary {(s,y):y = 0}, and the subscript on \i refers to that value 
of y. 

We calculate from P| 

— (t,0;s ,2/o) = - / <f* ~2^r [(* - s) + - 2/) ] L =3/=0 ^o(s) 

= 0- (8) 

As always in the boundary-integral method, we must make a careful distinction between 
the value of such an integral exactly on the boundary (x = 0) and the limit of the integral 
as the boundary is approached from the interior. For the latter we have by (J3J) and (jSJ) 
(97 



^(t,x;s ,y ) = - J ds—[(t-s) +x] fi (s) 



^-T^oOO asx|0, (9) 
and it is this object that must be chosen to satisfy the Neumann boundary condition: 

implies 



dG, . , dGo 9 7 

= — (t,0 + -s ,y =^ + tt 
ox ox ox 



and hence 



/io(t) = 2 ^(t, + ; s , y ) = ^[(f - s ) 2 + Z/o 2 ]^ (10) 



7(t,x;s ,yo) = -|^ / dsln[(t-s) 2 + x 2 ][(s-s ) 2 + l/oT 1 - <U) 



But from the method of images we know a more elementary formula for 7: 

7(t, x; s , y ) = G (t, x; s , -y ) = ~ -- ln[(* - s ) 2 + (x + y ) 2 }. (12) 
The equivalence of fill I) and (j!2j) is shown in Appendix A| 
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4. Case 2: Infinite star graphs 

Let d v be the number of edges meeting at the central vertex. On edge j there is a 
coordinate x, sometimes written Xj , equal to at the vertex. Without loss of generality 
we can take yo to be located on edge j = 1. 

In analogy with (El) and (J7|) we construct the Green function in the form 

/oo />oo 
dsGo{t,x;s,0)fij{s) + / ds — — (t, x; s, 0)^(s) 
oo j —oo 

l Z" 00 

= -— / gLs lnf^-sf + x 2 ]^.?) 

+ -/ dslit-sf + x 2 }- 1 ^). (14) 

(Here the dependence of [i and v on the vertex t> is suppressed, along with that on 

The Kirchhoff boundary conditions are 

Gj(t,0; s ,y ) = G k (t,0; s ,y ) for all j,k, (15) 
dv dG 

J2^(t,0;s ,y )=0. (16) 
j'=i 

They translate into 

7 j (t,0 + ;so,?/o) = 7fc(^,0 + ; s ,?/o) for all j,k ^1, (17) 
7i(t,0 + ; s ,yo) = 7fc(^,0 + ; s ,y ) - G (t,0; s ,y ) (for any k ^ 1), (18) 

The symmetries of these equations suggest the further ansatze 

All fij are equal, (20) 
All Uj are equal except v\ , (21) 

d v 

E "i = °- ( 22 ) 
j'=i 

Then (fTTj) is satisfied, and (fT8l) becomes, by analogy with ([9]), 

vi{t) - v h {t) = -2G {t, 0; s , y ) = ^- ln[(* - s ) 2 + y 2 }. (23) 
Sum ([23]) over all k ^ 1 and use ([22]): 

" 2 ln[(t - s ) 2 + yo 2 ] = (4 - l)^i ~ /^k = - y^^fc = dvVi ■ 
71 fc^i i=i 
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(l- j)^ln[(t-ao) 2 + W) a ], (24) 



and then 



= - 77—7- ln[(t - s ) 2 + y 2 ]- (25) 



Finally, to impose (1191) and find fij we note from (I14p that 



d'y- x f°° f°° 

-^■(t,x;s ,y ) = - — y ds [(t - s) 2 + x 2 ]~Vj(s) + y ds 



<9 2 G 



-00 

and hence from (TT9D and (122 



dxdy 



(t,x; s,0)vj(s) 



dG x f°° 

- — — (t, 0; s , y ) = - lim — / ds [(t - s) 2 + x 2 ]' 1 /^(s). 
dx xio 2tt y_ 00 ^ 

It follows by (JH]), ([ID, and §D that 

/^(t) = ^-[(t-«o) a + W?r 1 . (26) 

Inserting (|24l) - (|26j) into (fl4l) we arrive at 

1 p , ln[(t - s) 2 + x 2 } 



4:7r 2 d v J.oo (s - s ) 2 + Vo 2 

+ {*» - 5 J 4^ L ds x (t-.) 2 + * 2 • (27) 

The two integrals in (1271) are identical except for the interchange (t,x) (so>Z/o)- 
Moreover, at the end of the previous section, and in ( 1A.3I) . we observed that 

y f°° , ln[(t - s) 2 + x 2 } , r/ , 2 , .,. , , 

this function is actually invariant under that interchange. (In (T28l) it is assumed that x 
and ?/o are positive.) So we finally arrive at 

1 



7j (t, x; s , y ) = [Sji ~ f ) ^ H(t - s ) 2 + (x + y ) 2 ], (29) 

where we now allow for yo to be located on any edge, I. This is equivalent to known 
results for quantum star graphs (e.g., [19], Sec. 3B], [32], Ch. 3]). From ([3]), the formula 
for the cylinder kernel T is 

t-(^v)^ + ^t|t^ • (30) 

which is equation (35) of [9] with the Robin parameter a = 0. 
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5. Case 3: Y = (0, L) 

I write d/dn for an inward normal derivative in the usual sense of bounded domains. 
In the graph context such a derivative is outward from a vertex. 

In the present model the boundary has two parts, {(s, y): y — 0} and {(s,y):y = L}, 
the boundary condition is still Neumann, and the obvious ansatz is ([6]) with 



7(t,x;s ,2/o) = / ds G (t, x; s, 0)// (s) + / dsG (t,x;s,L)/j, L (s). (31) 

J — oo 

Evaluate the derivative: 



d^y 1 f°° x 

Tx {t, x; so, Vo) = J ds {t _ s)2 + x2 M 



1 f°° x — L 

-Y,j_J a {t-sY + ( x -LY ^ s) - (32) 

The limits of the normal derivatives from inside are 

|1(0+) S gft 0+; *, = , oW + JL £ ds _£_ (33 ) 

•97 _ d-f 1 1 Z" 00 L 

dn = ~ Yx^ L ~ ] So ' yo) = ~2 + 2^ L ds (t - a )2 + V ^ S) - (M) 

(The sign constraint in ([9]) becomes L—x J. in going from (|32l) to (l34l) .) The appropriate 
boundary conditions are 

) = -^(^°; s o,2/o), ) = +-=—(*, L;s ,2/ ). (35) 
on ox on ox 

From these relations we obtain the basic boundary-integral equations of the problem, 
lt\ 1 y " 1 1"^ , L 

. . 1 L - 3/o 1 Z" 00 , L 

»dt) = - (t _ So)2 + (L _ ^ + - y^d. (i _ s)2 + L2 /*(.). (36) 

More abstractly, the system (I3"6T) is 

w) = ~ n \7~r~Tf \? + - / rfs 7: m> ; r9 /M s )> ( 3T ) 



vr (t - So ) 2 + <*(t>, Z/o) 2 vr (f - sY + L 2 

where d(v, yo) is the distance from vertex v to yo > and U is the other vertex. Still more 
abstractly, fix a vertex v* ; then 

dG f dG 

/*«„(*) = 2 -Q^-{t,v*;s ,yo) +2 J _ -^-{t, v*; s, v)/i v (s), (38) 

where a sum over v is implicit in the integration over the total boundary dT = R x dT 
of the graph, and it happens, as in (|Sj), that 
9Gq 

— — (t, v*; s,f*) = 0. (39) 
on 

Even more abstractly, the integral equation has the form 

H = g + Kfx, (40) 
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so that formally 

H = (1 - K)- l g Q ~ (1 + K + K 2 + ■ ■ -)g . (41) 



Postponing to Appendix B the issue of the convergence of the Neumann series (141 p . 
we examine its zeroth term, which produces the first-order (not leading) term in the 
series for G. That is, one drops the integrals in (J36l) and substitutes (J36l) into (jSJ) and 
(EU): 

G(t,x;s ,y ) - G (t,x; s 0l y ) 

= - -L {ln[(t - So ) 2 + (x + 2/ ) 2 ] + ln[(t - s ) 2 + (2L - x - y ) 2 ]} , (42) 

where ( )28|) has been used (with L — x > 0, L — y > 0). Clearly, (j4"2"|) comprises the 
single-reflection terms in the standard solution by the method of images (images at — yo 
and 2L — y ) . 

The first-order term in (|4l|) should therefore yield the two-reflection terms in the 
image solution. Substituting the zeroth-order \x (the first terms in (|36|) ) into the integrals 
in (!36j) . one gets 



TT (t - s ) 2 + Vi ^ [(t - s) 2 + L 2 ) [(s - s ) 2 + (L — y ) 2 ] 

1 2L-y 

tt (f - s ) 2 + (2L - yo) 2 ' 

1 L-y 1 r° Lyo 

M*) " 77 72 TTTTTt 77X2 = 777 ds 



vr (t - so) 2 + {L- y ) 2 vr 2 J _ [(t - S ) 2 + L 2 } [(s - s ) 2 + y 2 } 

L + y 



(43) 



vr (t-so) 2 + (L + yo) 2 
(For the evaluation of these integrals see ( 1A.1I) .) The resulting additional terms in G 
are 

1 f r , , ru ^2 , 21 2L - y 



|y rfs ln[(t - s) 2 + X 



4tt 2 1 J.oo (s - s ) 2 + (2L - y ) 



2 



(44) 



+ £ rfshl|(t - s)2 + (i - l)2, (7^tT^} 

= - i- {ln[(t - So ) 2 + (x + 2L - y ) 2 ] + ln[(t - s ) 2 + (2L - x + y ) 2 ] } • 
Exactly as expected, these terms describe images at yo — 2L and yo + 2L. 



6. Case 4: General compact Kirchhoff quantum graphs 



Notation: v is a vertex of degree d v ; e is an edge. Whenever considering a fixed vertex v, 
we can assume that the edges are parametrized so that x = x e = at v; each incident 
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edge (e 6 £„) will then have a terminal vertex t e at which x e = L e . Let e$ be the edge 
containing the source point y Q . 

Ansatz (generalizing ( Tl~4l) and (EH)): 

G e = 5 eeo G + 7e . (45) 

For every v and every e G E v there will be charge and dipole densities fi e v{t) and u ev (t) 
(with hidden dependence on (s , yo) as usual). For temporary notational purposes, given 
e choose one of its two vertices, v, as its initial vertex. Then 

/oo poo 
dsG (t,x;s, 0)fi ev {s) + / dsG (t,x;s,L e )fi et (s) 
-oo J — oo 

/OO QQ r-OO QQ 

ds -^-(t,x;s,0)u ev (s) - J ds-^-(t,x;s,L e )u et (s). (46) 

Boundary conditions: Also for temporary purposes, consider a fixed v and assume 
the edges in E v are numbered 1, . . ., d v . (Of course, this list may or may not include 
eo •) Then the boundary conditions at v are, for all e and / in E v , 

7 e (*, 0; s , yo) +S eeo G (t,0;s ,yo) = 7/(t, 0; s , y ) + 5 feo G (t, 0; s , y ), (47) 

e=l e=l 

For (HHP we calculate 



^(t,x; s ,y ) = j ds^p-(t,x-,s,0)/j, ev (s)+ f ds^p-(t,x; s, L e )/i et (s) 

OX J —oo J —oo 

ds — ^° (t, x; s, 0)u ev (s) - J ds ^^ (t,x;s,L e )u et (s). (49) 

Integral equation system: Of the three simplifying symmetry relations (]20|) — fT22j) . the 
first and third are fundamental (essentially characterizing the solutions for the Dirichlet 
and Neumann subspaces of the graph — see, e.g., [2H E32] )■ They can immediately be 
generalized to the present situation: 

[lev = Hfv for all e, / e E v , (50) 

d v 

X>e, =0. (51) 

e=l 

However, (121]) was a symmetry that may no longer hold, because even if yo does not fall 
on any e E E v , some of those edges are "closer" to yo than others. As before, (1501) and 
( IBTj) assure that all nontrivial integrals at v in (j4"6"]) - (j4"5|) cancel. They also reduce the 
number of unknowns from 2d v to d v , which is the number of independent equations in 
( 147]) and (j48|) . For the trivial integrals we have as usual 

^-(t,x;s,0) = -^-[(t-s) 2 + (x-0) 2 }- 1 — > — | — s) asxjO. (52) 
Thus (T4T|) reduces to 

-v ev (t)+ I dsGo(t,0;s,L e )n et (s) - ds -^{t, 0; s, L e )u et {s) + 5 eeo G {t, 0; s , yo) 

^ J -oo J -oo oy 

= same with e — > f. (53) 
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This set of d v — 1 independent equations, together with the constraint (15"TI) . determines 
v ev (t) for all e 6 , if the v et and ft et are known. With the aid of (J50j) . fj4Bl) reduces to 
a single equation to determine /i et) (t): 



4 

2 



""^ POO QQ POO Q 2 Q 



dG 

dx 



(t,0; s ,y )^25 eeo 



(54) 



(The last factor in ( |54l) is equal to either or 1, unless eo is a loop at t>, in which case 
it equals 2.) 

So as in section [5] we have in (|53|) - (|54|) an integral equation of the form 

/2 = g + Kjl, (55) 

where \x is a vector- valued function whose d„ components are, for each v, any d v — l 
independent choices of the quantities v ev (constrained by (j5T]l ) and any one of the (equal) 
quantities /i e „ . Of course, in any concrete case the "et" notation needs to be resolved 
in favor of a fixed labeling of the vertices. In principle all integrals in the Neumann 
series can be formally evaluated recursively, and the results for Tit, x\ s, y) must be the 
same as those obtained by Wilson [32, 14J. Convergence of the series is not immediately 



obvious (see Appendix B). However, convergence of the series for the total energy has 
been proved in [6]. 
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Appendix A. Evaluation of integrals 

Lemma: Let x and yo be positive. Then 

ds 



oo 



[(t-s) 2 + x 2 ][(s-s ) 2 + y 2 ] 
(t — s) ds 



71 



x + Vo 



xy (t 

71 



so) + \ x 
t- s 



yo) 2 ' 



^ \(t - s) 2 + x 2 ] [(s - so) 2 + yo 2 ] y (t - s ) 2 + (x + y Q )< 



(A.l) 
(A.2) 
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Proof: Tediously, these integrals can be evaluated by standard methods — either 
residues or partial fractions. (The point (t, x) = (s , yo) requires special attention with 
a continuity argument.) 

Proposition: Let x and yo be positive. Then 

ln[(t 



ds 



- s) 2 + x 2 } 
so) 2 + yo 2 



71 



\n[(t - s y + (x + y ) % 



(A.3) 



[s - s Q y + y<f y 

Proof: Note that differentiating (1A.3I) with respect to x yields ( lA.lj) . while 
differentiating (1A.3[) with respect to t yields (1A.2H . It follows that the integral in (1A.3I) 
equals the right-hand side of (1A.3I) plus a constant, C, that is independent of x and t. 
The substitution s = s — s reduces the integral to 



00 d ~ s _ s ° ~ 5 ) 2 + x2 ] 

-00 



(A.4) 

from which it is easy to see that it is symmetric in t and So . Therefore, C is independent 
of so too. To fix C it suffices to consider the special cases x — > 0, t = 0, so = 0: 

In s 2 



C 



ds 



s 2 + y 2 



71 1 2 
— hi y 

yo 



1 

yo 



— / da 



2 -\ 


/•oo 


In a 


yo j 


—00 


a 2 + l 


4 


/ 


lncr 


yo 


10 


a 2 + l 










hi(aV; 

a 2 + 1 



do - 



da 



— lny 



In a 



a 2 + 1 



cr 2 + 1 



(A.5) 



because the substitution a — > 1/a converts the second term to the negative of the first. 
Remark: (IA.3I) is an instance of a formula given by Balian and Bloch [31 (7.2)], 

dG (ra) 



Go(rf 



da a 



dn r 



G (ar') 



(A.6) 



where r' is the image of r' in the plane S. They conclude it just by noting that in this case 
the boundary-integral/multiple-reflection expansion must agree with the image solution, 
assumed previously known. One of the main motivations of the present exercise was 
to verify flA.6j) directly, or, to put it differently, to derive the image solution from the 
(more general and fundamental) boundary-integral solution. The symmetry of fT25]) in 
(t, x) and (s, yo) is an instance of a symmetry of the right-hand side of (1A.6I) noted by 
Hansson and Jaffe ' Ji3 • (A. 12)]. 



Appendix B. Convergence issues 

Appendix B.l. Case 3 

Extrapolating from (H3|) . one can see that the series (14T!) will not converge, because the 
terms fall off only linearly with the distance to the image charge. (The corresponding 
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series for the Dirichlet case is conditionally convergent because of sign alternation.) The 
series for G, starting with (j42|) and (jHJ), is even worse: the terms grow logarithmically. 

In studying vacuum energy we are usually interested in derivatives of T. Their series 
converge better, because the differentiations build up powers of (t — s) 2 + (yo + ■ • -) 2 in 
the denominators. One can say that the original series converges distributionally with 
respect to test functions possessing sufficiently many antiderivatives that are also test 
functions (in the simplest case, test functions orthogonal to the constant functions). 

The standard test for convergence of a Neumann series like (|4ip is that the norm 
of K as an operator from some Banach space into itself be less than 1. In our case (see 
(I36l)-(l38l) the kernel function is built from 



dsK (t,s) = l = dtK (t,s). (B.2) 



so that 



It follows that Kq has norm 1 as an operator in either -^ 1 (R) or L°°(R). Furthermore, 
the iterated kernels have the same general form and the same norm; for example, 

/°° 1 2L 
ds K (t, s)K (s, So ) = - T+ , 2M , on2 (B.3) 
7T (t - s ) 2 + (2L) 2 

by (lA.lj) . So, in the multiple-reflection expansion we are operating right on the circle of 
convergence in the plane of a formal parameter multiplying K ( "marginal convergence" ) . 
Actually, one apparently needs flB.3j) to complete the abstract proof that the operator 
L 1 -norm is 1, not possibly something smaller. The foregoing remarks about conditional 
and derivative convergence show, however, that the latter possibility does not hold here. 
Note in passing that K is certainly not trace class or Hilbert-Schmidt, since 
ds Kq(s, s) is divergent. (The full operator K of (1371) has the matrix-valued kernel 




K = 

which, being totally off-diagonal, has a vanishing trace in the naive sense that 

ds K vv (s, s) = 0, 



E 



but that does not make it trace class.) The Fredholm solution is worse than divergent, 
because even the individual terms in the numerator and denominator series do not 
exist. Of course, this is to be expected for a convolution integral equation on an infinite 
interval. (This issue does not arise in the standard literature [H, [2], [3j [22J, [30J [31j [H [28] . 
which concentrates on compact boundaries.) If the graph were to be studied at finite 



temperature (Appendix C), the interval would be finite and the Fredholm theory would 



be applicable, possibly providing better convergence than that of the Neumann series. 
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Appendix B.2. Case 4 

The foregoing argument does not readily extend to the integral equation system ( l55l) . 
because the undifferentiated Go appearing twice in ( 1531) is not L 1 (nor L°°) in the "time" 
variables. 

At first glance it is not obvious even that the individual terms in the iterative 
solution will converge. Inspection of (J53H53I), together with (0}, ([A3]) , and f[Q|) . 
shows by induction, however, that each term in v ev {t) is logarithmic like Go and each 
term in /i e „(t) has behavior like dGo/dx as t — > oo. (There are always enough spatial 
derivatives to avert a disaster.) It seems likely that, as in Case 3, one is sitting on 
the edge of a region of convergence, and the convergence of 7 e can be improved by 
taking spatial derivatives. (See also the plausibility argument and numerical verification 
in [12] •) However, it is not possible to prove convergence by inspection without control 
over the number of terms in each order of the series and the magnitudes of their 
numerical coefficients. 

On the question of convergence of the series, some improvement is achieved by 
changing variables from v to v' . Equations fl53l) and fl54|) can be rewritten 

ds-^(t,0 ]S ,L e )fi et (s) 

/oo d 2 G dG 
ds s,L e )v et (s) + 2-7^(t,0; s ,y )5 eeo 

= same with e — > f, (B.4) 

r> <*« poo r, d v poo P)2ry 

^> ~d~ v ^J ~dx^ ' S ' Le >^ s > + d~ v ^J dx~dy~^ ' S ' Le)Vet ^ s > 

= ^-^(^,0;s ,i/o)^ 5ee o- ( B - 5 ) 

u e=l 

Observe now that 

&r<V ' 2tt [(t - s) 2 + (x - y) 2 } 2 " ds dt 1 J 

and 

dt 1 ' X ' *' yJ 2tt (i - s) 2 + {x- y) 2 ds ■ [ ' 

Some integrations by parts then put the system into the form 

/OO QQ 
ds ~0f (t,0;s,L e )fi et (s) 

/OO QQ QQ 

= same with e — > f, (B.8) 
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dv POO O s~y n dv 



2 \ f°° QQ 2 . \ /*°° QQ 

9 r)G ~ dv , 

~ JO (t,0;s ,y )J2See - (B.9) 



cL dx 



Recall that e and / are edges attached to the same vertex v, "t" refers to the other 
vertex of the edge in question, e is the edge containing the source point y , and the 
derivatives are in the direction away from v. Because of constraints fl50|) and ([5T]) . 
effectively there are Yl v d v independent components of v' and p, (twice the number of 
edges in the graph), although superficially the number of unknowns is twice that, and 
the number of equations even larger if all pairs {e, /} are counted. Again one can 
conclude by induction and Appendix A that all terms in the v' and p, series are of the 
same type as their respective source terms in (1B.8|) - (1B.9I) . 



Besides the redundancy of components, there are two technical obstacles to showing 
that (1B.8j) - (lB.9l) manifests an integral operator of norm at most 1. The first is that 



"time" derivatives of Go [( IB. 71) ] do not fall off as fast in the "time" variables as space 
derivatives [(jl])] do. However, one can hope to apply the generalized Young inequality 
[8], Theorem 0.10] with the weight function p(s) = (1 + |s|) _1 in all the integrals. With 
respect to that measure, every integral kernel in (1B.8j) - (lB.9j) is, with respect to the 



variables t and s separately, an L 1 function with norm at most 1 (in view of (IB. 21) and 
some trivial estimates). It follows that each term defines a mapping of norm at most 1 
in the Banach space LP (1 < p < oo) of functions such that 

|/(s)| p p(s) 1 ~ p ^ < oo (B.10) 

(if p = oo, the space of functions such that \ f(s)\/p(s) is bounded almost everywhere). 

The second obstacle is that each equation contains two such terms, not one, and 
after redundant components are eliminated from ( IB. 81) via ( 15T|) the number of terms 



will be still larger. Therefore, the obvious bound on the one-sided L 1 norms of the 
matrix- valued integral kernel is larger than 1, and one cannot conclude by this method 
that the series even marginally converges. 



Appendix C. Finite temperature 

The equilibrium state at temperature T should be described by replacing the usual 
cylinder kernel by a function periodic in t with period (3 = 1/T. In general, this 
construction can be implemented by replacing G by 

oo 

G T (t,x;s,y)= ^ G (t, x; s - N/T, y), (C.l) 

iV=-oo 

or by solving the partial differential equation with periodic boundary condition directly, 
or by summing the appropriate eigenfunction expansion with a Boltzmann weighting 
factor. 
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In total dimension 2, where G is given by (pQ), the literal (10.11) fails because of 
divergence. However, by making the arbitrary constants in ([1]) dependent on N one can 
arrive at 

Gt — — -— ln[cosh(27rT(a; -y))- cos(2nT(t - s))) + C, (C.2) 

47T 

whose first derivatives are 

dG T T sin(27rT(t- s)) 



<9t 2 cosh(27rT(x - y)) - cos(2vrT(t - s)) 

dG T T sinh(27rT(x - y)) 



(C.3) 
(C.4) 



dx 2 cosh(27rT(x - y)) - cos(27rT(t - s)) ' 

One can check the correctness of (1C.2j) by showing that it satisfies all the necessary 
conditions: 

(i) Periodicity in t with minimal period 1/T. 

(ii) Delta-function initial singularity: When t — s and x — y are small, Gt ~ Go , which 
is known to have the correct behavior. 

(iii) Asymptotic behavior at spatial infinity: logarithmic growth, just like Go • 

(iv) Partial differential equation 
d 2 G T + &Gr = Q 



dt 2 dx 2 
Moreover, as in (1B.6I) there holds 

g 2 G a 2 Go 
dx dy ds dt ' 
Remark: Alternatively, to deduce (IC.4I) from fIC.lj) . use the identity 

n2 + (h+ n\2 = ^-{coth[7r(a + ib)\ + coth[7r(a - ib)}} 



(C.6) 



a 2 + (b + N) 2 2a 

N=—oo v ' 



7T sinh(27ra) 



(C.7) 

a cosh(27ra) — cos(27r6) 

(The first version of this formula was listed in [10]. The simpler final form follows by 
elementary identities.) Going backward (a Mittag-Lefner expansion), one can get the 
analogous formula (1C.3I) . and then ( 1C.2I) is easy to guess. 
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